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Abstract 

Free  damped  vibrations  of  linear  elastic  structures  composed  of  uniform 
beam  elements  with  a  continuous  distribution  of  mass  are  studied.  Axial, 
torsional  and  flexural  vibrations  are  considered.  The  amount  of  damping, 
which  can  be  either  internal  or  external  viscous  type,  varies  among  the 
various  beam  elements  of  the  structure  resulting  in  many  critical  damping 
possibilities.  A  general  method  is  developed  which,  with  the  aid  of  dynamic 
stiffness  influence  coefficients  defined  for  every  element,  determines  the 
"critical  damping  surfaces"  of  the  system.  These  surfaces  represent  the 
loci  of  combinations  of  amounts  of  damping  leading  to  critically  damped 
motion  and  thus  separating  regions  of  partial  or  complete  underdamping  from 
those  of  overdamping.  The  dimension  of  a  critical  damping  surface  Is  equal 
to  the  number  of  independent  amounts  of  damping  present  in  the  system, 
while  the  number  of  these  surfaces  is  infinite,  i.e.,  equal  to  the  number 
of  degrees  of  freedom  of  the  system.  Three  examples  are  presented  In  detail 
to  illustrate  the  proposed  method  for  determining  critical  damping  and 
demonstrate  its  importance. 
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1.  INTRODUCTION 


The  importance  of  damping  as  a  means  of  reducing  the  response  of  a 
vibrating  structural  system  is  well  known.  Conventionally,  the  amount  of 
damping  in  a  linear  structural  system  is  expressed  as  a  percentage  of  the 
critical  damping  or  modal  critical  damping  values  depending  on  whether 
damping  is  everywhere  the  same  in  the  structure  or  varies  modally,  respect¬ 
ively.  Thus,  it  is  possible  to  estimate  directly  the  amount  of  damping 
in  the  structure  and  to  characterize  that  structure  as  underdamped,  over¬ 
damped  or  critically  damped.  This  knowledge  consequently  helps  one  to 
control  the  response  by  appropriate  changes  of  the  damping  in  the  structure. 

However,  for  linear  structural  systems  with  different  viscous  damping 
values  for  some  or  all  of  their  members,  the  problem  of  determining  critical 
damping  becomes  much  more  difficult,  because  many  critical  damping  possibili¬ 
ties  arise.  This  problem  is  of  considerable  importance  because  the  availa¬ 
bility  of  different  damping  values  for  different  members  of  a  structure 
provides  a  more  rational  way  of  representing  damping  properties  and  permits 
more  effective  response  control  by  taking  advantage  of  the  freedom  of  varying 
the  damping  of  a  large  number  of  elements. 

The  problem  of  critical  damping  is  part  of  the  general  problem  of 
structural  free  damped  vibration,  which  Is  concerned  with  the  determination 
of  natural  frequencies  and  modal  shapes  of  viscously  damped  linear  systems. 
Necessary  and  sufficient  conditions  under  which  discrete  and  continuous 
damped  linear  dynamic  systems  possess  classical  normal  modes  have  been 
established  by  Caughey  and  O'Kelly  [1].  In  a  recent  paper,  Beskos  and  Boley 
[2],  studied  free  viscously  damped  vibrations  of  linear  discrete  systems  in 
which  the  amount  of  damping  varied  among  the  various  structural  members, 
thus  resulting  in  many  "critical  damping  surfp-ses."  These  surfaces  represent 


the  loci  of  combinations  of  amounts  of  damping  leading  to  critically  damped 
motion,  and  thus  separate  regions  of  partial  or  complete  underdamping  from 
those  of  overdamping.  A  general  method  for  the  determination  of  critical 
damping  surfaces  of  linear  discrete  systems  was  developed  in  [2].  That  method 
is  extended  in  this  paper  to  certain  continuous  linear  dynamic  structural 
systems.  These  include  one,  two  or  three  dimensional  structures  consisting 
of  uniform  beam  elements  with  a  continuous  distribution  of  mass,  undergoing 
flexural,  torsional  or  axial  free  motion,  with  either  internal  viscoelastic 
or  external  viscous  damping.  The  method  developed  in  [2]  is  applied  to 
these  systems  in  conjunction  with  the  use  of  a  new  kind  of  dynamic  stiffness 
influence  coefficients  defined  for  the  aforementioned  motions  (flexural, 
torsional  and  axial)  on  the  basis  of  the  exact  solution  of  the  equation  of 
free  damped  motion  of  a  beam  element.  Thus,  the  dynamic  problem  is  reduced 
to  a  static-like  one  and  the  exact  solution  of  the  problem  is  obtained.  The 
use  and  importance  of  dynamic  stiffness  influence  coefficients  in  treating 
free  and  forced  vibration  problems  of  beam  structures  has  been  demonstrated 
elsewhere  [3-8]. 

To  the  authors  knowledge,  there  is  only  one  work  in  the  literature, 
namely  that  of  KolouSek  [9],  which  deals  with  viscously  damped  frameworks 
with  a  continuous  distribution  of  mass  and  different  amounts  of  damping 
among  the  various  structural  members.  However,  that  work  deals  only  with 
the  underdamped  steady -state  forced  vibration  case  by  employing  a  kind  of 
dynamic  stiffness  Influence  coefficients  In  convex  number  form.  Although 
only  a  certain  class  of  continuous  structures  is  considered  in  the  present 
paper,  namely  that  of  beam  structures,  the  results  obtained  are  representa¬ 
tive  in  that  they  demonstrate  special  features  contnon  to  all  continuous 
structures  characterized  by  an  infinite  number  of  degrees  of  freedom.  Other 
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continuous  structures  for  which  dynamic  stiffness  influence  coefficients 
can  be  constructed  can  also  be  studied  by  the  proposed  method.  If  this  Is 
not  feasible,  a  finite  element  discretization  and  modeling  of  the  structure 
as  a  discrete  system  with  a  finite  number  of  degrees  of  freedom  can  be 
always  done  and  the  method  of  [2]  then  applied.  Three  examples  dealing 
with  axial,  torsional  and  flexural  vibrations  are  presented  in  detail  in  this 
paper  to  illustrate  the  proposed  method  and  demonstrate  the  importance  of 

'  i 

critical  damping  surfaces. 

2.  FREE  DAMPED  VIBRATIONS  OF  A  BEAM  ELEMENT 

In  this  section,  dynamic  stiffness  influence  coefficients  for  free 
axial,  torsional  and  flexural  vibratory  motions  of  a  damped  linear  elastic 
uniform  beam  element  are  defined  and  constructed.  Either  internal  visco¬ 
elastic  or  external  viscous  damping  is  assumed.  Internal  viscoelastic 
damping  is  accounted  for  by  assuming,  for  reasons  of  simplicity,  that  the 
beam  material  is  a  Kelvin  solid,  i.e.,  with  a  one-dimensional  constitutive 
equation  of  the  form 

o  3  H  (1  +  g  d/dt)  c  ,  (1) 

where  a  is  the  stress,  c  Is  the  strain,  H  stands  for  the  modulus  of  elasticity 
E  or  the  shear  modulus  S,  g  Is  the  damping  coefficient  and  t  represents  time. 
More  general  viscoelastic  models  as  described,  for  example,  in  [10]  could 
also  have  been  used  without  any  particular  difficulty.  Equation  (1)  indicates 
that,  under  one  dimensional  states  of  stress,  the  formulation  of  the  damped 
beam  problem  can  be  obtained  from  that  of  the  corresponding  undamped  one 
by  simply  replacing  H  by  H(1  +  g  d/dt).  When  external  viscous  damping  is 
present.  It  is  accounted  for  in  the  displacement  equation  of  the  beam  motion 
by  a  damping  force  per  unit  length  P  proportional  to  the  velocity  and  opposing 
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the  motion,  i .e. , 

P  =  -  f(dv/dt)  ,  (2) 

where  f  is  the  coefficient  of  damping  and  v  is  the  beam  displacement. 

Consider  a  uniform  linear  elastic  beam  element  1-2  of  length  L  (Fig.  1) 
undergoing  axial,  torsional  and  flexural  free  damped  motions  which,  on  the 
basis  of  a  strength-of-materials  theory,  are  governed,  respectively,  by 
the  following  uncoupled  equations: 

EAu"  +  g  EAu"  -  mu  =  0  , 

(AC/J)4>"  +  g(AC/J)$"  -  m$  -  0  ,  (3) 

EIv'"'  +  g  EIv""  +  mv  =  0  , 

for  internal  viscoelastic  damping  of  the  Kelvin  type,  and 
EAu"  -  f u  -  mli  =  0  , 

(AC/J)<T  -  f$  -  m$  =  0  ,  (4) 

EIv""  +  fv  +  mv  »  0  , 

for  external  viscous  damping,  where  u  «  u(x,t),  v  *  v(x,t)  and  <|>  =  <J>(x,t)  are 
the  axial,  lateral  and  angular  displacements  of  the  beam,  respectively,  E  is 
the  modulus  of  elasticity  g  and  f  are  the  coefficients  of  internal  visco¬ 
elastic  and  external  viscous  damping,  respectively,  m  is  the  mass  per  unit 
length  of  the  beam.  A, I,  J  and  C  are  the  area,  moment  of  inertia,  polar 
moment  of  inertia  and  torsional  rigidity  of  the  cross-section  of  the  beam, 
respectively,  primes  Indicate  differentiation  with  respect  to  the  distance 
x  along  the  length  of  the  beam  and  dots  indicate  differentiation  with  respect 
to  the  time  t. 

By  assuming  solutions  of  the  form 
u  *  I  eXt  , 


<j>  *  <J>  eXt  , 
v  =  v  eXt  , 


(5) 


where  u,  <J>  and  v  are  functions  of  x  only,  and  X  is,  in  general,  a  complex 
number,  equations  (3)  and  (4}  are  reduced  to 
u”  -  K*u  *  0  , 

-  Kp  -  0  ,  (6) 

v'"‘  +  4K4v  «  0  , 

where 

K2  *  mX2/EA(l  +  gX)  , 

K2  =  mX2/(AC/J)(l  +  gX)  ,  (7) 

4K4  =  mX2/EI(l  +  gX)  , 
for  internal  viscoelastic  damping,  and 
K2  =  (mX2  +  fX)/AE  , 

K2  *  (mX2  +  fX)/(AC/0)  ,  (8) 

4K4  =  (mX2  +  fX)/EI  , 
for  external  viscous  damping. 

The  general  case  of  free  damped  motion  is  described  by  the  solution  of 
equations  (6),  namely 

u  *  D.  eKlX  +  F.  e"Klx  , 

♦  -  D2  eK2x  +  F2  e‘K2x  ,  (9> 

v  *  eKx  (DjCosKx  +  FgSinKx)  +  e~Kx(D4cosKx  +  F^sinKx), 

where  Dp  Fp  D2,  F2,  Dp  Fp  and  F^  are  constants.  Consider  the  beam 
element  1-2  of  Fig.  2  with  nodes  1  and  2,  which  has  one,  one  and  two  degrees 
of  freedom  per  node  for  axial,  torsional  and  flexural  motions,  respectively. 
Figure  2  shows  the  positive  directions  of  the  nodal  displacements  and  the 
corresponding  nodal  forces  for  the  three  kinds  of  motion.  The  dynamic 
stiffness  Influence  coefficient  Dij  Is  defined  as  the  force  at  the  1th 
degree  of  freedom  due  to  a  displacement  l.e**  at  the  j**1  degree  of  freedom, 
while  all  the  other  displacements  are  zero.  On  the  basis  of  this  definition 
and  by  using  the  displacement  functions  (9)  as  well  as  the  force-displacement 
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relations 

U(x)  =  AE  u' (x)  , 

T(x)*Ci',(x)  , 

V(x)  *-EIv"'(x)  , 

M(x)  =  -  EIv“(x)  , 

with  positive  directions  as  indicated  in  Fig.  1,  one  can  construct  the  Dij 
coefficients  for  the  three  kinds  of  motion  considered  here.  Thus,  with 
the  sign  convention  of  Fig.  2,  the  following  element  nodal  force-displace¬ 
ment  relations  in  terms  of  the  Dij  coefficients  result: 


for  the  axial  motion, 

,vi. k  Di2_j r 

hj  |_°21  022jl?2j 

for  the  torsional  motion,  and 


for  the  flexural  motion,  where 


Dil  =  °22  =  AEK1  coth^KiL)  * 

D-J 2  =  *  -  AEK1  cosech(K^L) , 


(12) 


(13) 


(14) 


°11  *  °22  *  CK2  coth(K2L)  * 

°12  “  °21  “  "  CK2  cosech(K2L^  » 


(15) 
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11 

8  °33 

*  2Nr  (y*  +  4Ysc  -  1)  , 

12 

8  °21 

■  -  °34  *  -  °43  '  NK^2  - 

2yO  -  2s2)+1  )  , 

13 

n 

o 

ii 

*  (4K2N/^)[-y2(s3  +  sc2  + 

c)  +  y(-s3  -  sc2  +  c2)]  , 

14 

8°41 

8  ‘  °23  8  "  °32  8  (2NK/W 

s3  +  sc2  +  s)(Y2  -  y)  , 

22 

8  °44 

*  N(y^  -  4ycs  -  1)  , 

24 

CM 

O 

II 

*  (2N/A)Cy2(s3  +  sc2  -  c) 

+  y(s3  +  sc2  +  c) 3  , 

N  *  2EIK/[y2  -  2y(1  +  2s2}  +  1]  , 

Y  *  e2KL  ,  s  *  sinKL  ,  c  =  cosKL  , 
and  where  ,  Kg  and  K  are  given  by  (7)  and  (8)  for  internal  and  external 
damping,  respectively. 


06) 


3.  CRITICAL  DAMPING  SURFACES 


Consider  a  three-dimensional  beam  structure  consisting  of  a  finite  number 
of  uniform  beam  elements  with  a  continuous  distribution  of  mass  and  with 
different  amounts  of  damping  under  free  motion.  The  damping  in  every  beam 
element  may  be  internal  viscoelastic  or  external  viscous,  or  a  combination 
of  these,  and  may  be, in  general,  different  for  different  kinds  of  free  motion. 
The  free  motion  of  every  element  (and  so  that  of  the  whole  structure)  is  a 
combination  of  axial,  torsional  and  flexural  motion.  Coupling  effects  between 
the  various  motions  are  neglected  in  this  work. 

The  general  equation  of  free  motion  for  the  above  described  beam  structure 
is  of  the  form  (e.g.  [7]) 

[D]{x}  =  [DH$}eAt  *  {0}  ,  (17) 

where  [D]  is  the  structural  dynamic  stiffness  matrix,  whose  elements  are 
combinations  of  the  dynamic  stiffness  influence  coefficients  Dij  for  axial, 
torsional  and  flexural  motion  given  by  (I4)-(16),  resulting  from  the  super¬ 
position  of  the  various  element  dynamic  stiffness  matrices,  and  {x}  and  {<!>} 
are  the  vectors  of  the  nodal  structural  displacements  and  displacement 
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amplitudes,  respectively.  The  matrix  [D]  is  thus  a  function  of  the  damping 
properties  of  the  various  beam  elements  and  the  eigenvalues  A  defined  by 
equations  (5).  Thus,  in  general, 

[0]  =  [D(hr,X)]  ,  (r  =  1,2, ...q)  ,  (18) 

where  h  represents  amounts  of  damping  and  q  is  an  integer,  in  general,  dif¬ 
ferent  from  the  order  p  of  the  matrix  [D]. 

Equation  (17)  has  nontrivial  solutions  if,  and  only  if, 

det[D(hr,A)]  =  |D(hr,A)|  =  0  .  (19) 

The  characteristic  equation  (19)  is  a  transcendental  equation  in  X  which 
has  an  infinite  number  of  roots  corresponding  to  the  infinite  number  of  degrees 
of  freedom  of  the  structure  under  consideration.  These  roots  can  be  negative 
real  or  complex,  resulting  in  overdamped  or  critically  damped  aperiodically 
decaying  free  motion,  or  in  underdamped  oscillatory  decaying  free  motion, 
respectively.  If  the  system  has  no  damping  (hp  *  0),  equation  (19)  is  satis¬ 
fied  for  an  infinite  number  of  imaginary  values  of  X  of  the  form  iu)Q,  where 
u>Q  represents  natural  frequencies  and  i  *  /T.  If  some  (all)  of  the  roots 
of  (19)  correspond  to  overdamping,  underdamping  or  critical  damping,  the 
structure  is  called  partially  (completely)  overdamped,  underdamped  or  critically 
damped,  respectively.  Once  the  roots  X  of  (19)  have  been  determined,  one  can 
solve  (17)  for  the  in  general  complex  modal  shapes  {$}. 


For  overdamping  or  critical  damping  the  roots  of  (19)  are  of  the  form 


X  =  -b  ,  (b>0)  , 

(20) 

and  (19)  becomes 

iO(hr,  -b)|  *0. 

(21) 

In  the  q-dimensional  space  with  coordinates  hr  (r  »  l,2,...q),  equation  (21) 
represents  a  family  of  q-dimensional  surfaces  corresponding  to  overdamping 
or  critical  damping.  The  problem  consists  of  determining  that  b  which  cor¬ 
responds  to  the  "critical  damping  surface."  There  are  actually  infinitely 
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many  "critical  damping  surfaces"  since  there  are  as  many  critical  damping 
possibilities  as  the  number  of  the  roots  A.  The  general  method  for  deter¬ 
mining  critical  damping  surfaces  of  linear  discrete  damped  systems  developed 
in  [2],  is  extented  here  to  linear  continuous  systems  described  by  (17). 

Thus,  having  in  mind  that  critical  damping  represents  the  threshold  between 

overdamping  and  underdamping,  one  can  conclude  that  among  the  S.  surfaces 

D 

described  by  (21),  the  critical  surface  Scr  is  the  one  for  which  the  damping 
is  a  minimum,  i .e. , 

(d/db)( |0(h  ,  -b)l)  =0  , 

r  (22) 

dh^db  *  0  ,  r  *  1,2,..,q  . 

In  principle,  equation  (22)  provides  the  bcr  as  a  function  of  the  hr's  and 
thus  the  equation  of  critical  damping  surfaces  is  given  by  (21)  with  b  =  bcr, 
i.e.,  by 

|D(hr,  -bcr)|  =  0  .  (23) 

In  practice,  however,  one  has  to  solve  the  system  of  simultaneous  nonlinear 
equations  (21)  and  (22)  numerically  in  a  manner  analogous  to  that  described 
in  [2]. 

The  above  method  for  determining  critical  damping  surfaces  is  quite  general 
and  applicable,  in  principle,  to  the  most  general  space  beam  structure;  how¬ 
ever,  the  practical  applicability  of  the  method  is  limited  to  small  order 
structures,  because,  to  the  authors  knowledge,  there  is  presently  no  efficient 
numerical  method  available  for  treating  equations  (22).  However,  things  are 
greatly  simplified  in  the  particular  cases  in  which  only  axial  or  torsional 
or  flexural  free  motion  arise  and  one  kind  of  damping  (internal  or  external) 
is  under  consideration.  Fortunately,  these  cases  deal  with  very  large  classes 
of  structures,  such  as  plane  frames  or  trusses  which  are  usually  considered 
to  undergo  only  flexural  free  vibrations  and  simple  and  composite  beams  or 
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shafts  undergoing  axial  or  torsional  free  vibrations,  respectively. 

Consider  a  beam  structure  undergoing  only  axial  or  torsional  or  flexural 
free  vibrations  under  conditions  of  either  internal  or  external  damping. 

The  coefficients  Klr’  K2r  and  Kr  for  axial,  torsional  and  flexural  free 
motion,  respectively,  of  the  rth  beam  element  (r  =  l,2,...,q,q<q)  are  given 
on  account  of  (7)  and  (8)  by 

=  [mr/(EA)r][A2/(l  +  grA)], 

V?2r  =  [mr/(AC/J)r][A2/(l  +  grA)],  (24) 

4Kj  =  [mr/(EI)r][A2/(l  +  grA)]  , 
for  internal  viscoelastic  damping,  and 

K^r  =  [mr/(EA)r](A2  +  28rA)  , 

4r  =  CVWJrlU2  +  26rX)  »  (25) 

4Kj  »  [mr/(EI)r](A2  +  28rA)  , 
for  external  viscous  damping  with 

fr  *  2mr8r’  (26) 

and  where  the  coefficients  gr  and  8p  in  (24)and(25)  are,  in  general,  different 
for  different  kinds  of  motion.  For  the  undamped  case  of  the  r**1  beam  element 
for  which 

Xn  =  1  “on  »  (n  ■  1*2,...,  •)  ,  (27) 

where  won  represents  the  natural  frequency  in  the  n^  mode  of  vibration  and 
i  *  /T  ,  equations  (7)  and  (8)  are  replaced  by 

KlrO  =  "  “on  ’ 

K2r0  ‘  ‘  “on  C"t/t*C/J)r3  .  (28) 

<o  '  -  “on  CV<E»r]  • 

When  one  particular  kind  of  free  motion  of  the  beam  structure  is  consid¬ 
ered,  the  structural  dynamic  matrix  [D]  consists  of  linear  combinations  of 
Dij  coefficients  with  the  same  type  of  K  coefficients.  In  that  case  (22) 


wpu'lfpuyp 
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becomes 

(d/db)( | 0 | )  =  Z(d/dK  )(|D|)-(dK/db)  =  0, 

r  r  (29) 

(dhr/db)  =  0  ,  r  =  l,2,...,q  , 

which  can  be  recognized  as  direct  extensions  of  the  corresponding  equations 
for  discrete  systems  [2],  The  derivatives  dK^db  for  the  various  cases  of 
motion  and  kinds  of  damping  can  be  computed  from  (24)  and  (25)  for  X  *  -b, 
and  in  conjunction  with  (29)2  take  the  following  forms: 
dKlr/db  =  Cmr/2Klr(EA)r  b(2-grb)/(l  -  grb)2  , 
dK2r/db  -  [m^K^AC/J)^  b(2-grb)/(l-grb)2  ,  (30) 

dK/db  =  [m/16^  (EI)r]  b(2  -  grb)/(l  -  gpb)2 


for  internal  viscoelastic  damping,  and 
dKlr/db  =  [mr/Klr(EA)r](b  -  0  r)  , 

dK2r/db  *  [mr/K2r(AC/J)r](b  -  8  p) ,  (31) 

dKr/db  -  [mr/8Kj  (EI)r](b  -  8  r)  , 
for  external  viscous  damping. 

The  particular  point  h^  =  h2=...=hq  =  h  of  a  critical  damping  surface 
for  which  conditions  (29)  hold  true  can  be  easily  determined.  Thus,  on 
account  of  (30)  and  (31)  with  8  *8  and  gp  »  g,  one  can  easily  see  that 

(29)1  is  satisfied  for 

2  -  gb  *  0  ,  (32) 

for  internal  viscoelastic  damping  and 


b  -  8  -  0  ,  (33) 

for  external  viscous  damping.  Substitution  of  (32)  and  (33)  into  (24)  and  (25) 
in  conjunction  with  X  =  -b,gp  =  g,  andBp  =8  and  comparison  of  the  results 
with  (28)  leads  to  a  critical  value  b  of  b  given  by 


b  =  oi 
ncr  on 


1,2,..., 


(34) 
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for  all  kinds  of  motions  and  damping  cases.  Thus,  (23)  takes  the  form 

lD(**r>  -  “0n>l  *°  (35) 

and  serves  to  determine  the  value  of  h  =  hn  for  every  n  numerically.  Notice 
that  in  this  case  the  structure  possesses  classical  normal  modes.  Equation 
(35)  also  indicates  that  a  beam  structure  with  a  continuous  distribution 
of  mass,  one  kind  of  damping  and  under  one  kind  of  free  motion  has  an  in¬ 
finite  number  of  critical  damping  surfaces  and  consequently  complete  critical 
damping  or  overdamping  are  achieved  for  n  -*■  «. 


The  above  results  could  have  been  obtained  in  a  different  way  as  follows 


When  h. 


*  hq  *  h,  combination  of  (24),  (25)  and  (28)  leads  to  the 


conclusion  that  the  value  of  A„  in  any  mode  n  will  be  such  that 

n 

A.2  /(I  +  gA  )  =  -oj2  (3 

n  s  n'  on 

for  internal  viscoelastic  damping  and 

A2  +  2  0A„  -  -  J  ,  (2 

n  n  on 

for  external  viscous  damping.  Under  conditions  of  underdamping 

\>  ■  -  bn  i  '“n  •  (: 

where  u  represents  the  damped  natural  frequency  In  the  0th  mode  and  by 
combining  (38)  with  (36)  and  (3Jj  one  obtains  the  relations 


(1/2)9w0„  , 


(l/2)mon(4  -  9 2<4,)'/2 


for  internal  viscoelastic  damping,  and 

bn’6'  “n  "  <"on  -  62)1/2  (40) 

for  external  viscous  damping.  Considering  critical  damping  as  the  threshold 

between  underdamping  and  overdamping,  one  can  Imagine  the  state  of  critical 

damping  as  the  limit  of  the  underdamping  state  as  the  amount  of  damping 

Increases  so  that  un  approaches  zero.  Thus,  the  conditions  of  critical 

damping,  in  view  of  (39)  and  (40),  are 


2^bncr  ’  bncr 


9 


(41) 


for  internal  viscoelastic  damping,  and 


6 


ncr 


b  =  cj 
ncr  on 


for  external  viscous  damping. 


(42) 

These  equations  are  identical  with  (32)-(34). 


Three  illustrative  examples  are  given  below. 

4.  EXAMPLES 


Example  1 

Consider  a  uniform  cantilever  beam  of  length  L,  mass  per  unit  length 
m  and  axial  rigidity  AE.  In  this  case,  damping  is  uniform;  the  example  is 
included  simply  to  illustrate  the  proposed  method  in  a  classical  case.  The 
equation  of  free  axial  vibration  of  this  beam,  obtained  from  (11)  for  = 

U2  *  0  and  TT|  =  0  reads 

°22  ^2  *  0  •  (43) 

Thus,  equation  (21),  on  account  of  (14),  (24)^,  (25)^  and  (43),  becomes 

coth(K.|L)  *  0  (44) 

where 

K2  =  (m/AE)[b2/(l  +  gb)3  (45) 

for  internal  viscoelastic  damping,  and 

K2  *  (m/AE)(b2  -  2Bb)  (46) 

for  external  viscous  damping.  The  equation  of  critical  damping  surfaces  is 
again  given  by  (44)  but  with  (45)  and  (46)  being  replaced  by 

K2  *  (m/AE)[u)2n/(l  +  9won)J  ,  (47) 

K2  =  (m/AE)(w2n  -  2  0uon)  ,  (48) 

for  Internal  and  external  damping,  respectively,  in  view  of  the  condition 
(34)  for  critical  damping.  The  natural  frequencies  won  can  be  obtained  from 
(44)  with  Kj  given  by  (28)^  and  are  of  the  form 

“on  *  (l/2L)(2n  -  1)  ir^ETm",  n  -  1,2,..,  (49) 

Thus,  for  both  kinds  of  damping,  the  critical  damping  surfaces  are  actually 
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an  infinite  set  of  points  which,  on  account  of  (32)  -  (34),  are  given  by 

9n  =  2/“on  (50) 

for  internal  viscoelastic  damping,  and 

=  ^on  (51) 

for  external  Viscous  damp!  no.. 

Example  2 

Consider  the  composite  shaft  1-2-3  of  Fig.  3  which  consists  of  two  uniform 
beams  of  circular  cross  sections  with  length  L  and  shear  modulus  G.  The 
polar  moment  of  inertia,  the  cross-sectional  area,  the  mass  per  unit  of  length 
and  the  external  viscous  damping  coefficient  are  symbolized  by  Jp  Ap  m^, 

Pj  and  Jjj,  Ajj,  mp.fip  for  the  beams  (1-2)  *  I  and  (2-3)  *  II,  respectively. 
The  equation  of  free  torsional  motion  of  this  shaft  is  obtained  with  the  help 
of  (12)  by  superimposing  the  dynamic  stiffnesses  of  the  two  beam  elements  I 
and  II  and  applying  the  boundary  conditions  ^  ^  =  0.  Thus,  equation  (21) 

becomes 

[(DJgJi  +  (DjVnl  =  0  (52) 

where,  from  (15)  for  a  circular  cross  section  (C  *  GJ), 

(°22)l  *  GJjK2Icoth(K2Il)  , 

(53) 

( Di* i ) 1 1  *  jCoth(K2j jt)  , 

with  K2j  and  K2I I  given  by  (27)2  with  X  *  -b  and  C  =  GJ.  Equation  (29), 
with  the  aid  of  (52),  (53)  and  (31)2  with  C  =  GJ,  yields 

(d(D22)l/dK2I){dK2I/db)  +  (d(°iVlI/dK2II)(dK2II/db)  =  0  »  (54) 

wKere  H  .2 

d  (1)22)1  ^Kzi  =  GiOjCcotKCKyL)-  KjjLcsck  (KaxL)l  > 

"HiiV*®1  stfeothfly.)  -  sch^K^n  , 

dK^/lb  =  (mi/K2iA1G)(b  -6i)  ,  1  -  I,  II  . 

The  natural  frequencies  of  the  system  can  be  obtained  from  (52)  with 


K2  ■  K2 
*21  *211 


'“on  <p's> 


» 


(56) 


-15- 


where 


P  5  (nij/Aj)  =  (nijj/Ajj)  ,  (57) 

represents  the  mass  density  of  the  shaft  material.  Use  of  (56)  and  (57) 
reduces  (52)  to 

(Jj  +  Jn)  cot(Kl)  *  0  (58) 

with  solution 


“on  =  0/2L)(2n-l)irv€7p,  n  =  1,2 . -.  (59) 

The  critical  damping  surfaces  are  described  by  (52)  with  b_„  obtained  as 
a  function  of  8  j  andBjj  from  (54)  and  (55)  and  form  an  infinite  set  of 
curves  in  theBj-Bjj  plane.  For  the  particular  case  of  Bjj  *  0,  Jj/Jjj  * 

2,  G/p  =  1  and  L  =  1,  a  computer  solution  of  the  simultaneous  equations  (52), 
(54)  and  (55)  provides  the  following  values  for  the  first  three  roots  Bj: 

Bj  -  2.03,  5.16,  8.31  .  (60) 


These  values  correspond  to  the  first  three  natural  frequencies 

“ol  *  */2  “  K571»  uo2  *  3ir/2  *  4'712*  “o3  *  5n/2  “  7'854  *  <61) 

and  indicate  that  the  critical  damping  surfaces  are  just  an  infinite  set  of 

points  along  the  Bj  axis. 

Example  3 

Consider  the  uniform  continuous  beam  1-2-3  of  Fig.  4  with  bending  rigidity 
El  and  mass  per  unit  of  length  m.  The  structure  consists  of  two  beam  elements 
of  the  same  length  L  but  with  different  external  viscous  damping  coefficients 
B  j  and  8g  as  shown  In  Fig.  4.  The  equation  of  free  flexural  motion  of  this 
continuous  beam  Is  obtained  with  the  aid  of  (13)  by  superimposing  the  dynamic 
stiffnesses  of  the  two  elements  (1-2)  *  I  and  (2-3)  ■  II  and  applying  the 
boundary  conditions  0^=  Vg*  vy*  0^  *  0.  Thus,  equation  (21)  becomes 

C ( ^44 ) j  +  ^22^11^  *  8  *  (62) 


I 
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where  (0^4)j  =  an<*  ^22^11  are  furnished  by  (lfi)  with  KI  and  Kn  9*v«n 

by  (25)^  with  X  *  -b. 

Equation  (29),  with  the  aid  of  (62),  (31 ) ^  and  (16),  yields 


(d(D22)I/dKI)(dKI/db)  +  (d^Jjj/dKjjHdKjj/db)  =  0 
dK./db  =  (m/8K?  EI)(b  -6.)  , 

d(D22)i/dKi  *  {[2(y2  "  4YCS  -  1)  +  8K.(y2  -  2ycs  -yc2  +  ys2)]  • 

•  [y2  -  2y(l+2s2)  +  1]  -  8K.(y2  -  4ycs  -  1 )(y2-y-2ys2-2ycs)>/ 

/  [y2  -  2y(l  +  2s2)  +  l]2  , 

Ki  =  (m/^  EI)(b2  -  2b  Bj)1''4  ,  i  =  I,  II. 

The  natural  frequencies  of  the  system  can  be  obtained  from  (62)  with 


(63) 


4KI  =  4KII  =  ’  (ra/EI)  “on  * 


(64) 


or  directly  from 


2D22  =  0  ,  (65) 

where  the  dynamic  stiffness  coefficient  without  damping  022  is  taken  from 
reference  [4]  or  [7]  and  is  of  the  form 

D22  *  EIK(sinKLcoshKL  -  cosKLsinhKL)/(l-cosKLcoshKL)  (66) 

with 

K2  =  u>  -vffiTTT.  (67) 

on 

Use  of  (66)  reduces  (65)  to 

tanKL  -  tanhKL  *  0  ,  (68) 

which  has  an  infinite  number  of  roots,  the  first  three  of  which  are  given 
to  third  decimal  accuracy  by 

KL  *  3.927,  7.069  ,  10.210.  (69) 

Thus,  by  combining  (67)  and  (69),  the  following  expressions  for  the  first 
three  natural  frequencies  are  obtained: 
u>on  •  (u^L2)/EUJr,  n  =  1,2,3, 

V2  -  15.421,  p2  =  49.971,  u2  *  104.244 


7 


I 


(70) 


The  critical  damping  surfaces  are  described  by  (62)  with  b  obtained 

as  a  function  of  8j  and  0jj  from  (63)  and  form  an  infinite  set  of  nonlinear 

curves  in  the  0j-3jj  plane.  All  these  curves  are  symmetric  about  the  line 

S  1=3 ii  ,  as  it  can  be  very  easily  seen  from  (16)  and  (25)j  with  X*-b  that 

an  interchange  of  0j  and  Sjj  leaves  (62)  unaffected;  this  was  expected 

in  view  of  the  symmetry  of  the  structure.  For  the  particular  points  of  the 

curves  for  which  8j*8jja8,  (42)  yields 

6  *  b  -  ui  .  (71) 

pn  ncr  "on  '  ' 

Figure  5  shows  the  first  three  critical  damping  curves  corresponding  to  the 
first  three  frequencies  of  (70).  These  curves  were  constructed  by  numerically 
solving  the  system  of  simultaneous  equations  (62)  and  (63)  for  El  »  L  *  m  *  1; 
the  computations  were  done  by  computer  using  complex  arithmetic  due  to  the 
fact  that  Kj  and  Kn  in  (63)4  are.  In  general,  complex  numbers.  It  was  ob¬ 
served  during  the  computation  that  the  values  of  bcr  for  every  curve  were 
very  close  to  the  value  of  bcr  corresponding  to  the  point 8 j  *Bjj,  i.e., 
to  the  value  8ncr  =  «Qn  .  This  suggested  an  approximate  construction  of  the 
curves  on  the  basis  of  (62)  with 
«?  ■  (m/EI)((uQn  -  2  6t»0„) 

(72) 

1  *  I,  II,  n  *  1,2,...,  °». 

These  approximate  curves  drawn  in  Fig.  5  for  El  *  l  »  m  *  1  practically  coincide 
with  the  "exact1'  curves,  Eventhough  this  suggests  a  simple  way  to  construct  very 
good  approximate  curves,  no  generalizations  beyond  the  above  observations  are 
presently  available.  The  critical  damping  curves  C^Cg, ..,(/»  of  the  structure 
of  Fig.  4  with  the  exception  of  the  C»  are  all  partially  critical  and  separate 


the  8  i  -  8  jj  plane  in  an  Infinite  number  of  regions  Rj,R£. ..,!*»  of  which 
only  the  first  four  are  shown  in  Fig.  5.  Region  R-j  represents  complete  under¬ 
damping,  while  complete  overdamping  Is  achieved  only  at  Infinity  as  all  the 
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other  regions  represent  partial  overdamping. 

The  damping  curves  of  Fig.  5  are  the  "exact"  damping  curves  of  the 
continuous  beam  as  computed  with  the  aid  of  the  dynamic  stiffness  influence 
coefficients  constructed  on  the  basis  of  the  exact  equations  of  motion. 

It  is  interesting  to  study  the  degree  of  approximation  obtained  by  con¬ 
structing  approximate  critical  damping  curves  of  the  continuous  beam,  on 
the  assumption  that  the  structure  is  discretized  by  lumping  the  mass  as 
shown  in  the  alternative  models  of  Fig.  6.  The  equation  of  free  flexural 
motion  of  the  two  degrees  of  freedom  model  of  Fig.  6(a),  on  the  assumption 
that  its  viscous  damping  matrix  is  diagonal  and  by  eliminating  rotational 
degrees  of  freedom,  takes  the  form 


where  x-j,  c^  and  x2>  cg  represent  vertical  deflections  and  damping  coefficients 
for  the  points  1  and  2,  respectively.  The  static  stiffness  coefficients  for 
a  beam  element  were  taken  from  [11].  The  natural  frequencies  u >Q  of  the  two 
degrees  of  freedom  system  are  the  roots  of  the  equation 


+  35 

35  -w*M’  +  115 

where 


(74) 


5  -  12EI/7U')3. 


(75) 


From  the  approximate  method  of  [2],  the  equations  of  the  critical  damping 
curves  of  the  system  are 


where 

C|  ■  2M'  6 1 »  c2  «  2M'  82,  n  »  1 ,2. 
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For  the  case  of  El  =  L  *  m  *  1, 


M*  =  0.5  ,  C  =  13.714  , 

(77) 

and  the  solution  of  (74)  is 

wQl  =  14.813,  wo2  =  19.596 

(78) 

forming  a  lower  bound  to  the  exact  solution  (70).  A  similar  computation 
for  the  four  degrees  of  freedom  model  of  Fig.  6(b)  provides  the  following 
values  for  the  natural  frequencies: 


wQl  *  15.357,  uo2  =  22.055,  u>o3  =  45.656,  wq4=  41.254.  (79) 

It  is  obvious  from  (72),  (78)  and  (79)  that  both  of  the  discrete  systems 
provide  an  acceptable  approximation  for  the  first  natural  frequency  only. 
Consequently,  both  of  the  discrete  systems  can  provide  a  good  approximation 
of  only  the  first  critical  damping  curve.  Figure  7  shows  the  first  critical 
damping  curve  for  the  continuous  beam  as  obtained  by  the  continuous  as  well 
as  the  two  discrete  models.  The  main  conclusion  of  this  comparison  study 
is  that  for  discrete  lumped-mass  models  only  the  first  few  critical  damping 
surfaces  are  close  to  the  exact  ones.  Accuracy  of  representation  of  the 
higher  damping  surfaces  Increases  by  increasing  the  order  of  the  discrete 
model . 

5.  CONCLUSIONS 

On  the  basis  of  the  preceding  discussion,  the  following  conclusions  can 
be  drawn  with  regard  to  the  free  motion  of  linear  structures  composed  of 
beams  with  a  continuous  distribution  of  mass  and  different  amounts  of  internal 
and/or  external  viscous  damping: 

1)  There  exist  infinitely  many  "critical  damping  surfaces"  for  every  struc¬ 
ture  which  represent  the  loci  of  combinations  of  damping  leading  to  partial 
or  complete  critical  damped  motion  and  thus  separating  regions  of  partial  or 
complete  underdamping  from  those  of  overdamping.  The  dimension  of  these 
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y  independent  amounts  of/ 

surfaces  is  equal  to  the  number  of  ^damping  present  in  the  system. 

2)  A  general  method  is  proposed  for  determining  the  equations  of  these 
critical  damping  surfaces  in  conjunction  with  dynamic  stiffness  influence 
coefficients  which  are  functions  of  the  damping  and  inertia  properties  of 
the  structure.  This  method  is  greatly  simplified  for  systems  which  undergo 
only  one  kind  of  motion,  namely  axial  or  torsional  or  flexural.  For  these 
systems  there  are  infinitely  many  partially  critical  damping  surfaces 
creating  one  region  of  complete  underdamping  and  infinitely  many  regions  of 
partial  overdamping  as  complete  critical  damping  or  overdamping  are  achieved 
only  for  infinite  damping  values. 

3)  The  critical  damping  surfaces  obtained  by  the  proposed  method  represent 
the  exact  solution  of  the  problem  since  they  are  constructed  with  the  aid 

of  dynamic  stiffness  influence  coefficients  based  on  the  exact  solution  of 
the  equation  of  motion.  For  complicated  beam  structures  or  other  structures 
with  continuous  distribution  of  stiffness,  mass  and  damping  properties,  such 
as  discs,  plates  and  shells,  a  discretization  of  the  structure  and  application 
of  the  methods  of  reference  [2]  is  probably  the  most  practical  approach.  The 
resulting  damping  surfaces  will  of  course  be  approximate,  with  a  degree  of 
accuracy  decreasing  as  surfaces  corresponding  to  higher  frequencies  are  con¬ 
sidered. 
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CAPTIONS  OF  THE  FIGURES 


Figure  1 
Figure  2 
Figure  3 
Figure  4 
Figure  5 

Figure  6 

Figure  7 


:  Positive  beam  displacements  and  forces  in  mechanics  convention. 

:  Positive  beam  nodal  displacements  and  forces. 

:  Composite  shaft  in  torsional  free  motion  for  Example  2. 

:  Continuous  beam  in  flexural  free  motion  for  Example  3. 

:  The  first  three  critical  damping  curves  of  the  continuous 
beam  of  Example  3. 

:  Lumped-mass  discretizations  of  the  continuous  beam  of 
Example  3. 

:  Approximations  of  the  first  critical  damping  curve  of  the 
structure  of  Example  3. 
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